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With bichromatic fields it is possible to deterministically 
produce entangled states of trapped ions. In this paper we 
present a unified analysis of this process for both weak and 
strong fields, for slow and fast gates. Simple expressions for 
the fidelity of creating maximally entangled states of two or 
an arbitrary number of ions under non-ideal conditions are 
derived and discussed. 

Pacs. 03.67.Lx, 03.65.Bz, 89.70,+c 



I. INTRODUCTION 

Quantum computing relies on the ability to perform 
a collection of unitary evolutions of a quantum system, 
composed of a number of two- level systems (the qubits) , 
and it is a key result that a small set of so-called univer- 
sal gates exists, which may form the basis for the entire 
computation 0] . The development of proposals for phys- 
ical implementation of quantum computing have followed 
different routes, according to the various views one may 
have on the quantum dynamical processes, (i): one may 
view a gate operation on a single or on several qubits 
as a controlled transition from the initial to the final 
states, and one may implement it by a Hamiltonian, or 
a sequence of Hamiltonians, that couple these states di- 
rectly, (ii): one may consider Hamiltonians that couple 
quite many states, but where unwanted operations are 
dynamically suppressed by resonance conditions or by 
'bang-bang' Hamiltonians |2| . (iii) : one may depart from 
a more systematic analysis of the Lie algebra generated 
(by commutation) from a given set of basic Hamiltonians; 
If one has access to Hamiltonians Hi and H2 with vari- 
able strength parameters K\ and k 2 , subsequent applica- 
tion over short time intervals dt of K\Hi, k 2 H2, —n\Hi 
and —K2H2, leads to the evolution operator (h — 1) 



^K 2 H 2 dt giKiHi dt ^—iK 2 H 2 dt^—iK.\H\dt 



3K1 K 2 [Hi,H 2 ]dt 2 +0 (dt 3 ), 



(1) 



so that effectively the Hamiltonian i[Hi, H2) is obtained. 
As expressed by Lloyd ||: 'By going forward and backing 
up a sufficiently small distance a large enough number 
of times, it is possible to parallel park in a space only e 
longer than the length of the car'. If Hi and H2 commute 
with the commutator [Hi , H2] , the higher order terms 
in dt vanish exactly and one may apply Hi and H 2 for 
arbitrarily large dt and 'make a round trip in the parking 
lot and park in one single operation'. 



The different proposals for quantum computing with 
trapped ions can be roughly categorized according to 
the lines above: In their original proposal Q|, Cirac and 
Zoller noted that lasers resonant with sideband excitation 
of the trapped ions couple the ground and first vibra- 
tional state conditioned on the internal state of the irra- 
diated ion, and subsequent irradiation of a second ion can 
couple its internal states conditioned on the vibrational 
state. We have formulated a proposal for two-bit || and 
multi-bit H gates in the ion trap, which makes use of 
resonance conditions to couple certain states of the two- 
particle system. In our proposal we apply bichromatic 
light which selects certain virtually excited intermediate 
states, and by choosing appropriate parameters we show 
that the desired internal state dynamics of the ions may 
be perfectly achieved, even if the vibrational degrees of 
freedom, used to couple the ions, are not in their ground 
state. Recently, Milburn Q has proposed a realization 
of a multi-bit quantum gate in the ion trap, which also 
operates when the ions are vibrationally excited: Ad- 
justing the phases of laser fields resonant with side band 
transitions, one may couple internal state operators to 
different quadrature components, e.g., position and mo- 
mentum, X and P, of the oscillatory motion. In Ref. Rj 
it is proposed to use the two Hamiltonians Hi — XiJ z P 
and H2 = \2JzX, expressed in terms of the collective 
spin operators = J2k^ k (£ = x iVi z \ where the sum 
is over the ions irradiated by the lasers, and where j^k is 
the spin operator for the atom k, which may be defined 
by the Pauli spin matrices j^k = &£k/2 (h = 1). By al- 
ternating application of the Hamiltonians Hi and H2 we 
may obtain the exact propagator 



e iH,T e iHir e - 



iH 2 T —WxT _ 



(2) 



because the commutator of the oscillator position and 
momentum is a number. The interaction contained in 
between the ions has been established via the vibra- 
tional degrees of freedom, but after the gate this motion 
is returned to the initial state and is not in any way en- 
tangled with the internal state dynamics. Milburn also 
considers the possibility of coupling different individual 
internal state operators successively to X and P, so that 
the commutator term provides the product of such oper- 
ators. 

In this paper, we shall demonstrate that our bichro- 
matic excitation scheme is in fact already a realization of 
the proposal by Milburn, and that gate operation more 
rapid than concluded in || is possible. We show that 
our bichromatic scheme implements a propagator of the 
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form c iA ( T ) J y which is analogous to the one obtained 
by Milburn @. In Ref. § it was shown that this propa- 
gator can be used to prepare maximally entangled states 
-^(\gg...g)+e l< ^\ee...e)) of any number of ions (N), where 
the fc'th letter denotes the internal state e or g of the 
k'th ion. These maximally entangled states, which have 
an interest in the their own right B, are produced by 
applying the unitary operator e ITr/ * to a string of ions 
initially in the state \gg...g), and they may be produced 
even without experimental access to individual ions in 
the trap. 

In this paper we focus on the preparation of maximally 
entangled states. This is both of convenience for the the- 
oretical presentation and to emphasize results which are 
most easily verified experimentally. However, the proce- 
dures described here also apply to quantum computation. 
With two ions illuminated by laser light the bichromatic 
scheme produces -^(\gg) — i\ee)) and together with sin- 
gle qubit rotation this evolution forms a universal set 
of gates which may be used to constuct a quantum com- 
puter. The control-not operation Jl]] for example may 
be obtained by applying single ion operations on each ion 
before and after the bichromatic pulse which creates the 
state -^{\gg) - i\ee)) from \gg). 

In section |l[ we recall our proposal for a two-qubit 
gate operation and we show that it is equivalent to the 
proposal of Milburn, with a harmonic rather than a stro- 
boscope application of Hamiltonian coupling terms. In 
experiments it may be difficult to fulfill the requirements 
for the analysis of sec. [H] to be precise, and in section 
III we address the fidelity with which certain entangled 
states may be engineered when we take into account the 
off-resonant couplings and the finite value of the Lamb- 
Dicke parameter. In sec. IV we study the influence of the 
environment on the system. We analyse the role of spec- 
tator vibrational modes and energy exchange between the 
ionic motion and thermal surroundings. A summary of 
our results and a conclusion are presented in section [y|. 



II. GATE OPERATION UNDER IDEAL 
CONDITIONS 

Ions in a linear trap interacting with a laser field of 
frequency lo may be described by the Hamiltonian 

H = H o + iJi n t 

H = v{a)a + 1/2) + uj eg ^ a zi /2 



(3) 



where v is the frequency of the vibration, a) and a are 
the ladder operators of the quantized oscillator, w eg is 
the energy difference between the internal states e and 
g, and Qi is the resonant Rabi frequency of the i'th ion 



in the laser field. The exponentials account for the po- 
sition dependence of the laserfield, and the recoil of the 
ions upon absorption of a photon. The positions of the 
ions Xi are replaced by ladder operators kxi = ^(a + a^), 
where the Lamb-Dicke parameter r\i represents the ra- 
tio between the ionic excursions within the vibrational 
ground state wavefunction and the wavelength of the ex- 
citing radiation. In Eq. (||) we have assumed that the 
laser is close to a sideband to « u> eg ± v for a single vibra- 
tional mode and that we may neglect the contribution 
from all other vibrational modes. We tune the lasers 
close to the center-of-mass vibrational mode where all 
ions participate equally in the vibration, so that the cou- 
pling of the recoil to the vibration is identical for all ions, 
i.e., Tji = T) for all i. For simplicity we also assume the 
same Rabi frequency for all ions participating in the gate 
fli = Q. In this section we will consider an ion trap op- 
erating in the Lamb-Dicke limit, i.e. the ions are cooled 
to a regime with vibrational numbers n ensuring that 
(ri + l)rj 2 << 1, so that we may perform the expansion 
e i v (a+ a i) ra 1+^(0 + a t). 



A. Weak field coupling 

In our previous article || we assumed that two ions 
in the string were both illuminated with two lasers of 
opposite detunings co — ui eg = ±5. With this choice of 
laser detunings the only energy conserving transitions are 
from \ggn) to \een) and from \gen) to \egn), where n is 
the quantum number for the relevant vibrational mode 
of the trap, cf. Fig. [p. We considered the weak field 
regime 77O << v — 6, where only a negligible population 
is transfered to the intermediate levels with vibrational 
quantum numbers n ± 1 . In this regime the effective Rabi 
frequency for the transition from \ggn) to |een) may 
be determined in second order perturbation theory 



(een\H int \m}(m\H in t\ggn) (nrj) 2 



6' 



(4) 



where we have used the intermediate states \m) = \egn±l) 
and \gen± 1), and where u m is the frequency of the laser 
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FIG. 1. Energy level diagram for two ions with quantized 
vibrational motion illuminated with bichromatic light. The 
only resonant transitions are from \ggn) to \een) (left) and 
from \egn) to \gen) (right). Various transition paths involving 
intermediate states with vibrational number n differing by 
unity are identified. 
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exciting the intermediate state |m). For the transition 
\egn) to \geri) we get the same effective Rabi frequency. 

The remarkable feature in Eq. (^) is that it contains no 
dependence on the vibrational quantum number n. This 
is due to interference between the paths shown in Fig. [l|. 
If we take a path where an intermediate state with vibra- 
tional quantum number n+ 1 is excited, we have a factor 
of n+ 1 appearing in the numerator (yn + 1 from raising 
and y/n + 1 from lowering the vibrational quantum num- 
ber). In paths involving the vibrational state n — 1 we 
obtain a factor of n. Due to the opposite detunings, the 
denominators have opposite signs and the n dependence 
disappears when the two terms are subtracted. The co- 
herent evolution of the internal atomic state is thus in- 
sensitive to the vibrational quantum numbers, and it may 
be observed with ions in any superposition or mixture of 
vibrational states. The coherent evolution may even be 
seen if the vibrational quantum number n changes during 
the gate due to heating Q . 



B. General field coupling 

We now consider the interaction without restricting the 
parameters to a regime where no population is transfered 
to states with different n. For this purpose it is conve- 
nient to change to the interaction picture with respect to 
Hq. In the Lamb-Dicke limit with lasers detuned by ±5 
the Hamiltonian becomes 

Hi n t = 2flJ x cos 5t — V2r]flJ y 

x [x(cos(^ — S)t + cos(f + 8)t) 
+p(sin(v- 5)t + sin(v + 5)t)], (5) 

where we have introduced the dimensionless position and 
momentum operators, x = -^(a+a^) andp = -^(a^—a), 

and the collective spin operators discussed above Eq. (||). 

Choosing not too strong laser intensities << S and 
tuning close to the sidebands v — <5 < < <5 we may neglect 
the J x term and the terms oscillating at frequency v + 5 
m Eq. (|), and our interaction is a special case of the 
Hamiltonian 

H int =f{t)J y x + g(t)J y p. (6) 

The exact propagator for the Hamiltonian (||) may be 
represented by the ansatz 

JJ(t) = e- iA ( t ) j2 y e ~ iF(t)JyX e~ iG( - t)JyP ', (7) 

and the Schrodinger equation i4zU(t) = Hi nt U(t) then 
leads to the expressions 



F(t) = / f(t')dt' 



G(t) 



A(t) 



g(t')dt' 
F(t')g(t')dt'. 



With f(t) 
— y/2rjri sin(i/ 

F(t) = 

G(t) = 
A(t) = 



— \Z2riil cos(v - 
5)t following from 

sin ( [v — d)t) 

v — o 



S)t and g(t) 
3) we get 



v-5 



[1 - cosO - 8)t)] 
1 



t - 



2{v - S) 



sin(2(> - S)t) 



(9) 



In the xp phase space the operator U performs transla- 
tions (x,p) — * (x + J y G(t) } p — J y F(t)) entangled with 
the internal state of the ions. 

Apart from a change of basis from J z to J y the in- 
teraction considered by Milburn may also be put in 
this form, with f{t) and g{t) alternating between zero 
and non- vanishing constants. Within the present for- 
mulation, the trick in Ref. is to use functions f(t) 
and g(t) such that F(t) and G(t) both vanish after a 
period r. At this instant the vibrational motion is re- 
turned to its original state and the propagator reduces 



to U(t) 



-iA{r)Jt 



we are left with an internal 



state evolution which is independent of the external vi- 
brational state. This decoupling is possible because the 
effective internal state transition is completed in the same 
amount of time for all vibrational components and be- 
cause the AC Stark shift of the atomic levels due to the 
laser fields are independent of the value of n. In the weak 
field case these properties are ensured by the interfering 
coupling amplitude in Fig. [I], see detailed discussion in 
Ref. Q. In the general case it follow from the formal 
structure of Eq. (Q). According to (||) the acquired fac- 
tor A(t) is equal to the area swept by the line segment 
between (G(i),0) and (G(t),-F(t)), as shown in Fig. |. 
If {G(t), —F(t)) forms a closed path, A(t) is plus (minus) 
the enclosed area if the path is traversed in the (counter) 
clockwise direction. In the proposal by Milburn succes- 
sive constant Hamiltonians proportional to x and p are 
applied and the area enclosed by (G(t),—F(t)) is rect- 
angular. In our proposal the area is a circle of radius 




(8) 



(G(t),-F(t)) 



FIG. 2. The paths traversed in phase space and the func- 
tion A(t) in Milburns proposal (rectangular) and in our pro- 
posal (circle). 
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\f2J y rflj [y — 6), as illustrated in Fig. 0. 

With the propagator in Eq. (0) we may calculate the 
time evolution of the system. Suppose that the ions are 
initially in the internal ground state and an incoherent 
mixture of vibrational state as described by the density 
matrix p tot = J2 n Pn\g--gn){g..gn\. The time evolution of 
the internal state density operator p = Tr„(/?'°') with any 
number of ions N may be found from p ai ... aN ,bi..b N (t) = 
J2n P n(g--g n \U 1i (t)\b l ..b N )(a l ..a N \U(t)\g..gn) (a,j,bj = e 
or g), where we have used J2 n \ n ) ( n \ — 1 to remove one of 
the summations over vibrational states. Here we list the 
density matrix elements for the case of two ions N = 2: 



Pgg,gg 



1 F(f) 2 +G(t) 2 



8 + 2° 



i,[ F{t)2+ 2 G{t)2 ) cos U(t) + \F{t)G{t) 



\ e -{F(t? + G(tr) Ln ^ F{t) 2 + G(t} 2^ 

3 1 F(t) 2 +G(t) 2 



E p > 

n 

xL 



Pgg,' 



SMI!) m (A, ) + l m o W ) 

+ i e -(nt) 2 +G(tf) Ln{2{F{t) 2 + G(t)2)) 



1(1 



-(F( t y+G(ty) 



xL n (2(F(t) 2 + G(t) 2 ))) 

i _ E11 f ± G 11 f T (F{t) 2 + G{t) 2 
"7T e 4 L n \ 



x sin ( A(t) + -F{t)G{t) 



(10) 



where L n is the n'th order Laguerre polynomium. 

These expressions can be evaluated in different 
regimes. In the weak field regime, ryfi << v — 6, the 
xp phase space trajectory is a very small circle, which is 
traversed several times. F(t) and G(t) are negligible for 
all times, and e - lF (t) J v x e - lG (t) J yP [ s approximately unity, 
such that we have an internal state preparation which is 
disentangled from the vibrational motion throughout the 
gate. Since A(t) w -ifVL 2 tj(y - 5) if {v - 8)t » 1 the 
time evolution corresponds to the one obtained from an 
effective Hamiltonian H = £IJ 2 , and Eq. ( JTo| ) describes 
simple Rabi oscillations between the states Jgg) and |ee). 
This is demonstrated in Fig. (a) which shows the time 
evolution described by Eq. (TO). The curves show sinu- 
soidal Rabi oscillation from \gg) to |ee) superimposed by 
small oscillations due to the weak entanglement with the 
vibrational motion. 

Outside the weak field regime the internal state is 
strongly entangled with the vibrational motion in the 
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FIG. 3. Time evolution of density matrix elements for two 
ions calculated from Eq (0). (a) Weak field regime (b) Fast 
gate. The first curve (counting from above at vt w 1000 
in (a) and vt w 130 in (b)) represents p 99l99 , the second is 
the imaginary part of p 99 , ee , the third is p ee ,ee, and the last 
curve is the real part of p 99 , ee . The ions are initially in the 
internal state \gg) and a thermal vibrational state with an 
average of 2 vibrational quanta. In (a) the physical parame- 
ters are S = 0.9u, r\ — 0.1, and fl — O.Iza In (b) the physical 
parameters are 8 — 0.95^, r\ — 0.1, and £1 — 0.177^. The pa- 
rameters in (b) are chosen such that a maximally entangled 
state -^{\gg} — *|ee)) is formed at the time vt w 250, where 
the circular path in Fig. ^ has been traversed twice. 



course of the gate. For successful gate operation we have 
to ensure that we return to the initial vibrational state 
at the end of the gate by choosing parameters such that 
G(t) = F(t) = 0, corresponding to (y — 8)t = K2ir, 
where K is an integer. A maximally entangled state is 
created if we adjust our parameters so that A(t) = —ir/2. 
This is achieved if the parameters are chosen in accor- 
dance with 



rift 



1 



2VK 



A" =1,2,3, 



(11) 



By choosing a low value of K such that an entangled 
state is created after a few rounds in phase space we may 
perform a faster gate than considered in the weak field 
case. See Fig. || (b), where we have used K — 2, and 
where a maximally entangled state -^(\gg) — i\ ee )) i s 
created at the time vt ~ 250. 

By combining the requirement ([ll]) with the condition 
(v — 5)t = K2tt we may express the time for the state 
preparation as 



<K. 



(12) 



In order to avoid off-resonant excitations of the ions we 
must require << 1 and rj 2 must be much less than 
unit y to f ulfill the L amb-Dicke approximation (see sub- 
sec. [II A and III B| ). For a given trap and/or laser in- 
tensity Eq. (|l2) sets a bound on the speed of the gate. 
In tabel | we give some numerical examples for the time 
of the gate for some typical experimental parameters. 
The CONTROL-NOT operation which plays a central role 
in quantum computation |Q may be created by a com- 
bination of single particle rotations and a bichromatic 
pulse with the duration described by Eq. p). The single 
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TABLE I. The time required to create the maximally en- 
tangled state -^(\gg...g) — i\ee..e)) with a Lamb-Dicke pa- 
rameter 77 — 0.1 for various trapping frequencies (v) and laser 
intensities (Q). The table shows the gate time if the entangled 
state is prepared after a single round in phase space. If the 
gate operation is accomplished after K rounds in phase space 
the time should be multiplied by \ // ~K. 
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particle operations may be performed much faster than 
the two qubit gates, so the time required to perform a 
control-not operation is also given by (|i~2]). 



III. NON IDEAL CONDITIONS 

In the previous section we used the Lamb-Dicke and 
the rotating wave approximations to arrive at an exactly 
solvable model. In this section we perform a more de- 
tailed analysis of the validity of the approximations and 
we estimate the effect of deviations from the ideal situa- 
tion in an actual experiment. The general procedure in 
the section, is to change to the interaction picture with 
respect to the simple Hamiltonian (JsJ) using the exact 
propagator in Eq. (|7]), and to treat the small deviations 
from the ideal situation by perturbation theory. The fig- 
ure of merit for the performance of the gate is taken to be 
the fidelity F of creation of the maximally entangled N- 
particle state \^ m ax) = l/v / 2(|g5"<?) — i\ee..e)), which in 
the ideal case is created at the time when A(t) = — tt/2, 
if the ions are initially in the \gg--g) state Q, i. e., 

in ax 

)■ (13) 



A. Direct coupling 

Going from Eq. (§) to Eq. (|) we neglected a term 
H d = 2flJ x cos(St). This term describes direct off res- 
onant coupling of g and e without changes in the vi- 
brational motion. For high laser power this term has a 
detrimental effect on the fidelity, which we calculate be- 
low. 

Changing to the interaction picture, we may find the 
propagator Ui(t) from the Dyson series 



Uiit) = 1 - i [ dt'H dJ (t') 
Jo 




dt'dt"H d: i(t')H dJ (t") + ... , (14) 



where the interaction Hamiltonian is given by H d j(t) = 
W (t)H d {t)U{t). Since H d {t) is oscillating at a much 
higher frequency than the propagator U(t), we may treat 
U (t) as a constant during the integration and we obtain 

oo 

UAt) = l-i— sm(dt)UUt) J x U(t) 
o 

~{l-cos(2St))U^(t)4U{t) + ... . (15) 

Near the endpoint, U(t) w e l (' K /' z ) J l and we obtain the 
fidelity 

F*l- -^(I-cob^t)), (16) 

where N is the number of ions participating in the gate. 
We plot in Fig. || the product of the fidelity due to 
the carrier (|l^) and the population of the EPR-state 
~^{\gg) ~ A ee )) expected from the time evolution in 
Eq. @. The result agrees well with the result of a 
numerical integration of the Schrodinger equation with 
the Hamiltonian (|5|). 
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FIG. 4. Population of the EPR-state A^Qgg) — i\ee}) near 
the optimum. The full line is obtained by a numerical in- 
tegration of the Hamiltonian Jcfc and the dashed line is the 
product of the expression in ( |lq ) and the expression for the 
fidelity obtained from Eq. ( |lo| ) . The parameters are the same 
as in Fig. | (b). 

If the duration of the laser pulses can be controlled 
very accurately in the experiment, so that one fulfills 
both ( fTl| ) and 2<5r = 2iG7r the effect of the direct cou- 
pling vanishes. If one cannot perform such an accurate 
control, the net effect of the direct coupling is to reduce 
the average fidelity by ^S- (=0.03 for the parameters of 
Fig. |). 
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B. Lamb-Dicke approximation 



In section |l| we used the Lamb-Dicke approximation 



-^(a-j-a^ 



1 + irj{a + a') to simplify our calculations. 



Now we investigate the validity of this approximation. 
In the weak field case, we can use the exact matrix 

elements (n\e tTl ^ a+a ^ \n + 1) = i?y S-^=j (r; 2 ), to obtain 

the effective Rabi frequency between \ggn) and |een) 



2 1 



n 



I - n 1 i2n+\) + rf[ h -n 2 + h -n + ~ 



(17) 



where f2 is given by Eq. (|]), and where L\ are the gen- 
eralized Laguerre polynomials 



m=0 



(-1)' 



(18) 



The effective Rabi frequency is no longer independent of 
the vibrational quantum number n, and the internal state 
becomes entangled with the vibrational motion, resulting 
in a non-ideal performance of the gate ||. 

To illustrate the effect of deviations from the Lamb- 
Dicke approximation, we consider again the production 
of an EPR-state A^{\gg) - i\ee)). With an n-dependent 
coupling strength the fidelity is 



F 



1 1 

n=0 



(19) 



where P n is the initialpopulation of the vibrational state 
n. We show in Fig. the evolution of the fidelity pre- 
dicted by Eq. ( JlQ| ) and obtained by a direct integration 
of the full Hamiltonian in Eq. (||) . Due to the deviation 
from the Lamb-Dicke approximation the effective Rabi 
frequency is reduced, c/., Eq. (|l7|), and the optimal gate 
performance is achieved with a duration that is longer 
than 7r/(2f2). The spreading of the values of Q n , causes 
entanglement with the vibrational motion which reduces 
the fidelity. With the parameters in Fig. |5| the maxi- 
mally obtainable fidelity is 0.92 obtained after a pulse of 
duration r i=s 1.9/0. 

With more than two ions, the time evolution of the 
system may be obtained by expanding the initial state 
\gg...g) on eigenstates of the J y operator: 



\gg--g) 



H) 



2 N/2 



N N 



fc=0 



\ 




\M y = N/2-k). (20) 



In the J y basis the propagator (Q) is diagonal and in the 
weak field regime (F(t), G(t) sa 0) with n dependent 
coupling strengths we get the fidelity 




FIG. 5. Evolution of the population of the EPR-state 
-j=(\gg) — i\ee}) for a vibrational thermal state with an aver- 
age of 5 vibrational quanta and r\ = 0.20. The dotted line is 
the prediction from Eq. ( ]l9| ) and the solid line is the result of a 
numerical integration of the Hamiltonian (|^) with parameters 
Q, — 0.02v and S = 0.9v. The discrepancy between the two 
curves at fit > 2 is due to additional off-resonant couplings 
which may be taken into account by multiplying the coupling 
strength by (dashed curve) gg. 



F 



n=0 




e i(N/2-k) 2 (ir/2-n n t) 



(21) 



In the limit of many ions (N » 1) and near the opti- 
mum (Q n t rs 7r/2) we may approximate this expression 
by assuming that k is a continuous variable and by replac- 
ing the binominal coefficient by a Gaussian distribution 
with the same width. In this limit the fidelity becomes 



F 



n=0 



1 



N(N-l)(Tr/2-n n t) 2 
4 



(22) 



1 



Expanding this expression to lowest order in rj and ad- 
justing the pulse duration to take into account the reduc- 
tion in the coupling strength we find to lowest order in 

V 



F = 1 - 



TT 2 N(N - 1) 



77 4 Var(n) 



at the optimum time 



(23) 



(24) 



where n and Var(n) are the mean and variance of the 
vibrational quantum number. 

In Eq. (^2|) and (^3|) we have replaced a quantity N 2 
following from the Gaussian approximation to (|2l| ) by 



G 



N(N-l). With this substitution (§|) and (§|) describe 
the fidelity well for all values of N. With the parameters 
of Fig. j| Eq. (||) yields F = 0.88 which is in good 
agreement with the numerical result in the figure. 

The equations ( p~T| - p3| ) were derived for weak fields, but 
they also provide an accurate description of the system 
outside this regime. To show this we note, that with 
Dichromatic light, H lnt in Eq. (p[) may be written as 



H, 



2f2 cos 



(St) J x cos (^]V2(x cos(vt) + psm(vt))^J 
(r]V2(x cos(vt) + psin(i4))) (25) 



- J y sin 



in the interaction picture with respect to Hq. An expan- 
sion of the trigonometric functions in this Hamiltonian 
leads to Eq. ([3|) which formed the basis of the discussion 
in section [n} The term proportional to J x is suppressed 
because it is far off resonance. The lowest order contribu- 
tion of this term was treated in the previous section, and 
we shall now consider corrections to the J y term which 
may have significant effects. In the interaction picture 
with respect to the lowest order Hamiltonian (^|), x and 
p are changed into x + J y G(t) and p — J y F(t) and to 
lowest non-vanishing order in rj the interaction picture 
Hamiltonian is 



H 3 =r 1 3 J y 



V2n 

12 



cos((v — S)t)hi(x,p) 
+ sm((v - 5)t)h 2 (x,p) , 



where 



hi (x, p) = 3a; 3 + xp 2 + pxp + p 2 x 
h,2(x,p) = 3p 3 + px 2 + xpx + x 2 p, 



(26) 



(27) 



and where we have used that F(t) and G(t) are propor- 
tional to rj. To calculate the effect of the Hamiltonian in 
(pq) we note that the propagator 



U: 



unt(*) = e[ 



-hl(x,p) 



h 2 (x,p) 



(28) 



is consistent with the Hamiltonian (|2^) until order rf , i.e. 
i dU %f {t) = (H 3 +0(rj 6 ))U 3 ,int- (But the full Hamiltonian 
contains terms of order rf and rf which are not taken into 
account in Uz; m t- These terms are included below). We 
are interested in the propagator at times r = K2ir / (v—8) 
where the vibrational motion is returned to the initial 
state. At these instants the exponents in Eq. ( p8| ) vanish 
and the propagator reduces to U 3 (t) — 1 such that it has 
no influence on the internal state preparation. 

Expanding the Hamiltonian to order rf we obtain the 
propagator to the same order in rj in the interaction pic- 
ture with respect to Ho in (|J) 

U e (r) = e -i^Jl[l-ri 2 {2n+l)+ V i (in 2 +in+^)] 



xc 



valid at times r = K2ir/ (v—5). The first exponential pro- 
vides the time evolution with the modified effective Rabi 
frequency in Eq. (|l7]). If we evaluate the propagator ( |29| ) 
in the weak field regime, the last two exponentials both 
vanish in the limit of large K when the req uirement ( |Tl| ) 
is inserted, and the time evolution in (29) is consistent 
with Eq. ( p7| - p^ ). The last two exponentials are also of 
minor importance for a different reason: In Eq. ([l7]) rf 
appears in the combination rj 2 n, whereas it appears as 
i] 2 in the last two exponentials of (||^) when the condi- 
tion ( |TT| ) is inserted. In situations where deviations from 
the Lamb-Dicke approximation are important rfn ~ I, 
the deviation is typically caused by a high value of n 
rather than a large value of rj (rj 2 << 1). In this case 
one may neglect the last two exponentials and the effect 
of the non-Lamb-Dicke terms are the same as in the case 
of weak fields as described by Eqs. ([l^ - ^3|) . To achieve 
the optimum operation of the gate with the parameters 
of Fig. H we have to ensure fir « 1.9 and there is a small 
correction to the condition in Eq. (O). 



V 2(v-<5) J 



(29) 



IV. EXTERNAL DISTURBANCES 

So far we have considered a system described by the 
Hamiltonian (^|), where only the center of mass motion 
is present in the ion trap and where the coupling of this 
mode to the surroundings is negliblc. In this section we 
shall remove these two assumptions and consider the de- 
crease in fidelity due to the presence of other modes in the 
trap and due to heating of the center of mass vibrational 
motion. 



A. Spectator vibrational modes 

With N ions in the trap, the motional state is de- 
scribed by 3./V non degenerate vibrational modes. With 
a proper laser geometry or if the transverse potential is 
much steeper than the longitudinal potential, the cou- 
pling of the laser to transverse modes will be neglible and 
the only contribution is from the N longitudinal modes. 
With iV vibrational modes the ion trap may be described 
by the Hamiltonian 

H = Hq + H ln t 

N 

H = vi(a\ai + 1/2) + uj eg ^ <j zl /2 

1 = 1 i 



where v\ and a\ and a; are the frequency and ladder op- 
erators of the i'th mode. The excursion of the i'th ion in 
the I'th mode is described by the Lamb-Dicke parameter 

rji^i which may be represented as 77^ = rj Nbi , where rj 



7 



and v refer to the center of mass mode as in the previous 
sections, and where b\ obeys the orthogonality conditions 

Eti Wi = and E^i b% = 0. ' 

The center of mass mode (I — 1), which is used to 

create the entangled states of the ions, has bj = 1/ >/~N 
for all ions and is well isolated from the remaning N — 1 
vibrational modes i^>i > ^/iv, so that we could neglect 
the contribution from the other modes in the previous 
sections. In this section we shall extimate the effect of 
the presence of the spectator modes. They have both 
a direct effect, due to the off resonant coupling to the 
other modes, and an indirect 'Debye- Waller' effect |Tl| l 
because the coupling strength of the center of mass mode 
is reduced due to the oscilations in the spectator modes. 
Below we shall calculate the direct and indirect effects 
separately. 

The lowest order contribution of the direct coupling 
to the spectator modes may be found by expanding the 
exponentials as in Eq. (|^). 



N 



2ttJ x cos St + > 6; [xifi (t) + Pl9l (t)] , (31) 



1=1 



where fi(t) = — \/2rjVl^ v / vi[cos{vi — 5)t + cos(i^ + S)t] 
and gi(t) = —^/2r]Q>^/v/vi[sm.{vi — S)t + sin(^ + 5)t], 
and where the internal and external state operators are 
defined by 9; = Eilo^W and Xl = 7a( a ' + a \) and 
pi = — a i)- Since the ladder operators for different 

modes commute, we may find the propagator for this 
Hamiltonian using the steps that lead to Eq. (Q) 



N 



U(t) =l[Ui(t), 



where 



Ui(t) 



e -iA l (t)&f e -iF l (t)e l x le -iGi(t)e lPl 



(32) 



(33) 



with the functions Fi, Gi and Ai defined analogously to 
Eq. (||). Note, that this is an exact solution of the Hamil- 
tonian ( |3l|) without the J x term, so that to lowest order 
in the Lamb-Dicke parameter it includes all effects of the 
coupling to the other modes. 

From the definition of @i it is seen that 9i = J y and 
the propagator U\ reduces to Eq. (fjj) in the rotating wave 
approximation. The other N — 1 propagators in ( [32] ) 
cause a reduction of the fidelity due to the excursion into 
the xipi phasespace of these modes. Expanding the expo- 
nentials, using (gg...g\&iQi> \gg...g) — SijiN/4 and S w v, 
and averaging over time we find 



o 2 N 



1=2 



vf/v 2 + 1 

(l/f/l^ - l) 2 



(34) 



where ni is the mean vibrational excitation of the I'th 
mode. 



In addition to the direct coupling to the spectator vi- 
brational mode, the fidelity is also reduced because the 
coupling strength is dependent on the vibration of the 
other modes. Unlike the direct coupling discussed above, 
this effect is not suppressed by the other modes being far- 
off- resonant, and it may have an effect comparable to the 
direct coupling. 

Due to the vibration of the ions the cou- 
pling of the i'th ion to the sideband is reduced 

from ir\\Jn + 1 to (7Jiri2...nN|e J Ei=i' )i >'( a!+a f -f 
ln 2 ...n N ) « i7 7 V^TT(f - + 1/2)). With 

this reduced coupling strength the effective propagator 
at times r = K2it / '(y — S) may be described by 



U(t) 



-iA(r)A 2 



where A = J27=i3vA l - EZiVtMi + 1/2)). In the 
Cirac-Zoller scheme || , the n-dependent AC Stark shifts 
caused by coupling to other vibrational modes lead to de- 
coherence, unless these modes are cooled to the ground 
state. In our bichromatic scheme, these internal state 
level shifts depend much less on the vibrational excita- 
tion. By expanding ( |35| ) around the optimum A(t) ps 7r/2 
we calculate the lowest order reduction in the fidelity 

p = 1 _ ^NjN-l) ^ y, Var(n,) 



(35) 



8 



N 



16 



i,l,V=l 



(6i) 2 (6f) 2 -l/iV i 



vivv/v 2 



nin v . (36) 



The expressions in Eqs. ( J34| ) and ( pq ) may be simpli- 
fied if the vibrational motion is in a thermal equilibrium 
at a given temperature. In a thermal state Var(n;) = 
nf + fii 1 nini' = nifiy for I 7^ I', and n\ < hiv/vi, and 
using these expressions we find the lower estimate for the 
fidelity 



02 

F > 1 - ifN— (ni<7i(JV) + a 2 (N)) 
for the direct coupling (|34|) and 

F > 1 _ ^ jV (^- 1 ) T? 4 ( -2 ff3(JV) + niai{N)) 
o 

n 2 (N-2) 2 



(37) 



16 



?? 4 (nfa 5 (^)+n 1 a 6 (iV)) (38) 



for the Debye- waller coupling (p6[) ; where the sums 

ai...ae may be derived from Eqs. ( p4|) and (|36"1). For 
example 03 (iV) = E;=i With the mode functions 

and frequencies of Ref. jn| these sums are readily eval- 
uated, and the results are shown in Fig. ^|. From the 
figure it is seen that 0-5,06 << 0-3,(74, so that the last 
term in Eq. ( |§8| ) may be neglected. All the sums have a 
very rapid convergence and we may estimate the fidelity 
by replacing the sums with their large N values, i.e. 



8 



tt 2 

F> 1 - T] 2 N— 0.8(7ii + 1) 
for the direct coupling (p3) and 



F > 1 



ir 2 N(N - 1) 



?7 4 (1.2n? + 1.4n x ) 



(39) 



(40) 



for the Debye- Waller coupling (pq). 

We note that Eq. © is derived from terms beyond the 
Lamb-Dicke expansion and it incorporates the reduction 
of fidelity due to deviations from the Lamb-Dicke ap- 
proximation in the center of mass mode, cf. the formal 
similarity of Eq. @ and Eq. (§|). 




10 20 30 40 50 60 70 80 90 100 

N 

FIG. 6. Evaluation of the sums ai...ae for different number 
of ions N. Starting from above at N ~ 5 the curves represent 

(74, C3, Cl, f2, (T6, and 1T5. 



B. Heating of the vibrational motion 

An ion trap cannot be perfectly isolated and the vi- 
bration of the ions will be subject to heating due to the 
interaction with the environment. Relaxation due to the 
interaction between the vibration and a thermal reservoir 
may be described by the master equation 



^P = -i[H,p] 



(41) 



where C(p) is of the Lindblad form 

£(P) = - \ E ( C l C ™P + pCLC m ) + £ C mP Cl (42) 



with relaxation operators C\ = \/r(l + nth)a and C2 = 
\/T(rith)a? , where V characaterizes the strength of the 
interaction, and nth is the mean vibrational number in 
thermal equilibrium. 



We calculate the effect of heating assuming that the 
ions remain in the Lamb-Dicke limit. Changing to the 
interaction picture with respect to the Hamiltonian (0), 
the time evolution of p is entirely due to the heating, 
i.e., the Lindblad terms which are transformed using the 
propagator (Q) 

G(t) - iF(t) 



Ci = tfdU = y/r(l + n th )[ a + J y 



C 2 =tfC 2 U= ^H h at + J, 



V2 

G(t) + iF(t) 
V2 



(43) 



The density matrix is most conveniently expressed in 
the basis of J y eigenstates, and by tracing over the vibra- 
tional states we find the time derivative of the internal 
state density matrix in the interaction picture 



J t PM y M- = -(M y - M' y fT{\ + 2n th 



G{t) 2 +F{tf 



■PMyMl 



(44) 



This equation is readily integrated, and at times r = 
K2n/ {y — S) we get 



PM y ,M'M = PM y ,M'(0)e-(^-^) 2m ^" 



(45) 



The initial state is expanded on the J y eigenstates as in 
Eq. ( pp| ) and the population of the initial state (which is 
ideally constant in the interaction picture) equals 



N N 



F 



2 2N EE 



j=0 k=0 




(46) 



For two ions this expressions can be readily evaluated 

(47) 



F= 3 - + 1 -e- 
8 2 



r (l+2"th) , 
4.K ' 



1 , r(i+2„ tfe ) 

-e 4K 
8 



In the limit of many ions (N >> 1) and short times 

( Tl ^ 1+ 4j( th ' > T <<; 1) we ma y agahi approximate the ex- 
pression in Eq. (Eq) by assuming that j and k are contin- 
uous variables and by replacing the binomial coefficients 
by Gaussian distributions with the same width. In this 
limit the fidelity becomes 



F = 



1 



1 



N r(l+2n th ) 



(48) 



For 2 ions the deviation between (|47|) and (|48|) is less 
than 0.02 for all values of F larger than 0.5. 

In the above expressions we have assumed the Lamb- 
Dicke approximation. This corresponds to a situation 
where the heating is counteracted for example by laser- 
cooling on some ions reserved for this purpose. If the ions 
are not cooled the heating will proceed towards high vi- 
brational numbers with a heating rate Trith and the heat- 
ing will eventually take the ions out of the Lamb-Dicke 
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limit. With strong fields [K ~ 1) the reduction in the fi- 
delity described by Eq. (^g) will ruin the entangled state 
before the heating has made a significant change to the 
vibrational state (Tn t hT > 1). For weak fields (K >> 1) 
however, the situation is different. With weak fields one 
may produce an entangled state even though the time 
required to entangle the ions is much longer than the de- 
coherence time of the vibrational motion which is used to 
communicate between the ions, i.e. if K > NTn t hT the 
effect of heating is small even though the change in the 
average vibrational number Yn t hT is larger than unity 
HD- Since the effective Rabi-frequency has a small de- 
pendence on the vibrational quantum number n as de- 
scribed in Eq. (|lj) the heating will have an indirect effect 
on the internal state preparation. This can be modelled 
by changing the probabilities in Eqs. (|l9|-p2|) into time 
dependent functions P n (t) reflecting the change in the 
vibrational motion occurring during the internal state 
preparation. 



or more quantum systems is a key ingredient of quantum 
information processing. The algebraic property (|2|) which 
allows coupling and temporary entanglement with such 
an ancilla may find wide applications in many different 
systems for quantum computation with different ancillac 
(photons, phonons, Cooper-pairs, etc.). However, oper- 
ators with a constant non-vanishing commutator (which 
allows the formal step from Eq. to Eq. (||)) only exist 
in infinite-dimensional Hilbert spaces Jl4| . In addition to 
the implementation in cavity QED realizations of quan- 
tum computing p^-p7[| where quantized cavity fields play 
the role of the vibrational modes, it thus seems very rel- 
evant to investigate to which extent the ideas underlying 
Eq. (Q) can be generalized to ancillae with a finite num- 
ber of states and, e.g., for communication across a linear 
qubit register by only nearest neighbour interaction. 
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V. CONCLUSION 

We have in this paper evaluated the possibility for 
preparation of entangled states of ions by illumination 
with bichromatic light. We have identified two regimes: 
(i) a weak field regime where single photon absorption is 
suppressed and where two-photon processes interfere in 
a way that makes the internal state dynamics insensitive 
to the vibrational state, and (ii) a strong field regime 
where the individual ions are coherently excited and the 
motional state is highly entangled with the internal state 
until all undesirable excitations are deterministically re- 
moved towards the end of the interaction. 

We have presented analytical estimates for the fidelity 
of the internal state preparation. These expressions are 
summarized in table [nj. The expressions for the fi- 
delity may be readily applied to experimental param- 
eters and they show that several ion trap experiments 
today are in a position to apply our proposal directly. 
In fact, using our proposal the NIST group at Boul- 
der has been able to produce the maximally entangled 
state -j^{\gggg) — i\eeee)) with four ions Q. In this 
experiment the heating of the center-of-mass mode was 
so strong that this mode could not be used to communi- 
cate between the ions. Instead the experiment used an 
asymmetric mode where all ions have the same amplitude 
but a different sign, i.e. \r]i\ are the same for all ions i. 
Apart from the center-of-mode such modes only exist in 
ion traps containing two or four ions, and the experiment 
could not go beyound four ions. In other existing traps 
the heating is much less significant [Q, and these traps 
may be employed to produce entangled states with more 
particles. 

The use of ancillary degrees of ff edom (center-of-mass 
position and momentum) to communicate between two 
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J x term in (I 



Nil 1 
~25 T 
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Lamb-Dicke vibration towards 

(n|e^ Q+at >| n+ 1) Spectator vibrational modes vibrational number 

n t h with rate Tnth 



^ irj^/n + 1 



(i) Direct coupling ,..s „ , „ T 
2 J .i j (n) Debye- Waller 

to other modes 



4 * 2 N(N-1) 
I 8 



4 7r 2 Af(Af-l) 
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